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Abstract 

The colour-suppressed tree amplitude in non-leptonic B decays is particularly sen- 
sitive to perturbative and non-perturbative corrections. We calculate the two-loop 
(NNLO) vertex corrections to the colour-suppressed and colour-allowed tree am- 
plitudes in QCD factorization. Our results are given completely analytically, in- 
cluding the full dependence on the charm quark mass. We then update theoretical 
predictions for a range of interesting observables derived from tttt, irp and pp final 
states that do not depend (significantly) on penguin contributions, and hence are 
now available with NNLO accuracy. We observe good agreement with experimental 
data within experimental and theoretical errors, except for observables involving 
the 7r°7T branching fraction. 
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1 Introduction 



Non-leptonic B decays are among the primary observables at the B factories. Their 
sheer number of more than a hundred final states offers a large number of observables: 
branching ratios, CP asymmetries, polarizations, as well as certain well-motivated com- 
binations thereof. The large dataset accumulated by the B factories and the Fermilab 
Tevatron therefore allows for a rich phenomenology and for precise determinations of 
CKM and unitarity triangle parameters. In the future LHCb and a possible SuperB fac- 
tory will further increase the amount of data; and they have the potential to discover yet 
unobserved channels and allow for an even more precise determination of quark flavour 
parameters. 

Despite the fact that non-leptonic B decays are theoretically not as clean as for in- 
stance rare or radiative B decays, the amount of available experimental data justifies the 
need for precise theoretical predictions. A successful attempt to gain control over the 
complicated QCD effects in non-leptonic B decays is QCD factorization (QCDF) [1-4]. 
Within this approach hadronic matrix elements are represented as convolutions of per- 
turbative objects (hard-scattering kernels) with non-perturbative quantities (light-cone 
distribution amplitudes, LCDA). It therefore constitutes a framework that systemat- 
ically disentangles perturbative from non-perturbative physics. The factorization for- 
mula, valid at leading order in Aq CD /mb and to all orders in a s , reads 



(M 1 M 2 \Q i \B) = im\ jif Ml (0) jf du T*(u) f M ^M 2 (u) 



+ / du dudv Z/ 7 (^,t>,w) f B 4>B(u) fM^Mxiv) fM 2 <pM 2 {u) >■ (1) 





The hard-scattering kernel T 1 starts from 0(1) and contains all vertex interactions, 
whereas T n starts from 0(a s ) and comprises hard spectator-interactions. On the basis of 
the flavour structure of the operators one distinguishes several "topological amplitudes" , 
which are known completely to next-to-leading order (NLO), i.e. through 0(a s ). The 
importance of accounting for radiative corrections to the hard-scattering kernels can be 
clearly seen from the expression for the colour-suppressed tree amplitude, which at NLO 
reads [1] 

a 2 (nn) = 0.220 - [0.179 + 0.077 z] NLO + [^] {[0.114] LOsp + [0.067] tw3 } . (2) 

One recognizes the large cancellation between the tree-level and the one-loop vertex 
correction (denoted by "NLO"), which is due to the fact that the Born term is small be- 
cause of colour suppression. This suppression is lifted at one-loop due to the appearance 
of additional colour structures^ The last term proportional to r sp is another colour- 
unsuppressed NLO correction that originates from hard spectator-scattering, which turns 



1 This suppression can be lifted only once, and all following terms should be natural in size as expected 
from a perturbative series. Hence, the size of the one-loop correction does not indicate a breakdown of 
perturbation theory. 
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out to be significant as well, and opposite in sign compared to the one-loop vertex cor- 
rection. 

From this discussion it is clear that the colour-suppressed tree amplitude is particu- 
larly sensitive to the next-to-next-to- leading order (NNLO), which is one of the reasons 
for going beyond NLO. This is also motivated by phenomenology since for certain quan- 
tities one observes deviations from QCDF predictions with experiment, especially in 
decay channels which are dominated by «2, such as the branching ratio and direct CP 
asymmetry of B — > 7i tt° [5]. It is therefore interesting and desirable to explore whether 
the NNLO QCDF corrections turn the amplitudes into the right direction to cure dis- 
crepancies between theory and experiment. Moreover, since direct CP asymmetries only 
start from 0(a s ), the 0(a 2 ) terms constitute only the first correction and are therefore 
needed in order to decrease the scale uncertainty stemming from the leading-order terms. 

Part of these 0(a 2 ) contributions is already complete, namely the one- loop hard 
spectator corrections to T , for which both the tree [6-8] and penguin [9, 10] topologies 
have been worked out. Later, the imaginary part of the vertex kernel T 1 [11, 12] was 
computed. There is clearly a demand for the complete 0(a 2 ) terms since leaving out 
one of the two terms in the formula (pQ) in a given order may lead to large unphysical 
effects, as is evident at NLO from (jSJ). In the present paper we compute the two- loop 
NNLO corrections to the topological tree amplitudes ot\ and «2, including the real and 
imaginary part, in a fully analytic form. Recently, the so far missing real part also 
appeared in [13]. What we add, besides confirming a technically highly non-trivial 
calculation, is an analytic expression for the charm and bottom mass dependence arising 
from massive quark-loop insertions into the gluon propagator. This should be useful for 
the NNLO calculation of the remaining penguin amplitudes that we plan to address in 
a later work [14]. We also discuss the numerical impact of the two-loop correction on 
B — > mi, 7rp, pp branching fractions and branching fraction ratios, updating some results 
of [5], which are now available with NNLO accuracy. A complementary analysis of these 
modes at NNLO can be found in [15]. 

Our article is organized as follows. In Section [2] we introduce the theoretical frame- 
work and explain the matching of QCD onto SCET. In Section [3] we describe the tech- 
niques applied during the two-loop calculation. In Section H] we give the master formulas 
for the hard-scattering kernels, and Section [5] contains the results of their convolution 
in Gegenbauer moments. The non-leptonic B — > tttt, Tip, pp branching fraction and ratio 
analysis is found in Section [61 We will conclude in Section [7J The lengthy analytic 
expressions for the unintegrated hard-scattering kernels and a discussion of the two-loop 
master integrals can be found in two appendices. 

2 Theoretical framework, matching QCD to SCET 

The starting point of our calculation is the effective weak Hamiltonian where the top 
quark and the heavy gauge bosons are integrated out [16]. In order to compute the 
0(a 2 ) corrections to a% and «2 we need the current-current operators of the Hamilto- 
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nian. Throughout this paper we work in the Chetyrkin-Misiak-Miinz (CMM) operator 
basis [17] since it allows to consistently use the naive dimensional regularization (NDR) 
scheme with anticommuting 75. We then have 

7U/f = ^= V: d V ub (dQi + C 2 Q 2 ) + h.c. , (3) 

where 

Q x = uYT a (l - l5 )b d lfl T a (l - l5 )u , 

Q 2 = uY(l-l 5 )b J Ta4 (1-7 5 )«. (4) 

Since we will use dimensional regularization the operator basis has to be extended by 
so-called evanescent operators. They are only non-vanishing for D 7^ 4 space-time di- 
mensions, and are needed in order to close the operator basis under renormalization. In 
our case they read 

E W = uYYYT a (l - i 5 )b d ltxlvlp T a {\ - l5 )u - 16Qi , 
E? = u<f>ff(l - l 5 )b d 7M 7,7p(l - 1> - 16Q 2 , 

E[ 2) = uYYYl a l X T a (l - 7s)& d^ liy%1 ^ x T a (l - lb )u - 20E; i} - 256Q 1 , 
E® = uYl u Yl a l X (l - 7s)& dw»7 P 7vlx(l - l 5 )u - 20E® - 256Q 2 . (5) 

Our computation amounts to performing a matching calculation from QCD onto 
SCET [18-21], which in the context of spectator scattering has been worked out in [6]. 
Adopting the notation of that paper for the calculation of the vertex correction means 
that we determine the matching coefficients of SCET four-quark operators [22]. There 
are two cases to consider depending on how the quark flavours from the operators (j4j) flow 
into the final state. The first possibility is that the d and u quarks from the (g?u)v_a 
current move in nearly the same direction, and are therefore described by the same 
type of collinear fields. This gives rise to the colour-allowed topological tree amplitude 
ai(MiM 2 ) , in which the second meson M 2 carries the flavour quantum numbers of [ud]. 
We refer to this as the "right insertion" of the operators (j4]) [6], and express the matrix 
elements of Qi as 

{Qi)=Y,Hia{O a ). (6) 

a 

The basis of non-local SCET four-quark operators on the right-hand side is given by 

7s)x ^+(l-7s)^, 

ishlilihix - 75)7±m 7 7±/m<A • (7) 



O x = X 



o 2 = x 



o 3 = x^(i 
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Our conventions are as follows. Meson Mi, which picks up the spectator antiquark 
from the B meson, moves into the direction of the light-like vector n_. The collinear 
quark field for this direction is £, with r/_£ = 0. The second meson M 2 moves into 
the opposite light-like direction n + . The collinear quark field for this direction is Xi 
satisfying rjl + x = 0. The effective heavy-quark field is labelled by the time-like vector 
v = (n+ + nJ)/2 with v 2 = 1 as usual, and satisfies tfh v = h v . The perpendicular 
component of any four-vector, including Dirac matrices, is defined by 

= n+ -p— + n_ -p-± +p^. (8) 

The SCET operators Oi are actually non-local on the light-cone, the n + -collinear part 
being of the form x(tn-)[. . -]x(0). The matching coefficients Hi a in are therefore 
functions of a single variable, usually chosen to be the Fourier conjugate of t, and the 
product H ia (O a ) must be interpreted as a convolution product. The non-locality, as well 
as the collinear Wilson lines that render the non-local product x(tn_)[. . -]%(0) gauge- 
invariant, and also the quark flavour, are not indicated explicitly in ©. We refer to [6] 
for further technical details concerning these issues. 

The second possibility, called "wrong insertion", arises when the u quark from the 
(m6)v-a current and u quark from the (Ju)v-a current move in nearly the same direction. 
The corresponding amplitude is the colour-suppressed tree amplitude 02(^1 M 2 ) and 
meson M2 is now made up of [uu]. To obtain this amplitude, the QCD operators are 
matched to SCET via 

(Qi) = Via (O a ) , (9) 

a 

where now the basis consists of 

Oi = Ct±(i - 7s)x x(i + l5h±aK , 

O3 = ^7^7±7l7l7l(l - 7s)x xi 1 + l^lXal^el^slx^l^K . (10) 

Note that the % fields, which carry the flavour of M 2 , now stand in two different fermion 
bilinears, hence the name "wrong insertion". The remarks above on non-locality apply 
to these operators as well. 

The SCET operators in ([7]) and (JTU1) are constructed such that all operators with 
indices 2 and 3 are evanescent, i.e. they vanish in D = 4 dimensions. Moreover, 0\ is 
Fierz equivalent to 0\ in D = 4 dimensions. Hence, we treat the difference 0\ — 0\ as 
another evanescent operator, and regard 0\ as the only physical SCET operator. Due 
to the absence of soft-gluon interactions between collinear fields in different directions 
in the leading-power SCET Lagrangian after a field redefinition [19], 0\ factorizes into 
a (xx) part, whose SCET matrix element is a light-cone distribution amplitude, and a 
(£h v ) part, which defines the soft part of a heavy-to-light form factor [23,24]. 
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Figure 1: Examples of factorizable (left) and non-factorzable (right) diagrams. By def- 
inition, in the former gluon lines do not connect the horizontal with the sloped lines 
(quarks that go into M2), in the latter they do connect them. The bold line denotes the 
bottom quark. 



On the QCD side of the matching relations ([6]) 
expansion of the renormalized matrix elements as 



we can write the perturbative 



(Qi) 



+ (S) 2 



4(1) , 7 (i) 4(0) , 7(1)4(0) 



A 



(2) , 7 (1) .(1) 7 (2) .(0) 7 (1) ,(1) 7 (2) .(0) 



+ *<2 Z^Af} + Z^A% + (-<) 5m^ A£A + (O a >W (11) 



for the right insertion, and for the wrong insertion we simply replace A — > A and 0^0. 
In (11 ip . and throughout the paper, a s denotes the strong coupling in the MS scheme 
with five active quark flavours at the scale /x, and a superscript on any quantity S will 
label the order in a s according to 



S 



S<(°) + «i s-d) + 

47T 



47T 



+ 0(a 3 ) 



(12) 



The renormalization factors Zy, Z Q , dm, and Z e2; t account for operator, coupling, mass, 
and wave-function renormalization, respectively. We will give more details on them in 
Section HJ The AfJ in (fTTj) denote bare £-loop on-shell matrix elements of the operators 
Qi in the effective weak Hamiltonian. For I = 2 we only need i — 1, 2 , but for £ < 2 we 
need in addition the tree-level and one-loop matrix elements of the evanescent operators 
due to terms such as 4- i?' in (jTTi) . The amplitudes A^p can be further split up into 
factorizable and non-factorizable diagrams, see Figure [T], according to 



A) 



(0 



A {e)f + A 

■"-in I Si 



(£)nf 



(13) 



which turns out to be convenient since in the matching procedure the factorizable dia- 
grams will cancel to a large extent. 

On the SCET side of the matching relation the matrix elements of O a assume a 
similar form to (jlip . but there is no mass counterterm for a HQET quark, and the 
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coupling is renormalized in the four-flavour theory (denoted by a hat on the coupling 
and renormalization constant), 



(Oa) 



M« + Y f 



(1) 



ext 



M. 



( 2 ) _j_ y(i)M (1) + Y {2) 



(2) a* + y 6 2 y« + z«m« 



ext °ab 



ab 



+ o(a2)} (o 6 ) 



(0) 



(14) 



The amplitudes are now labelled M^, and Y are the SCET renormalization constants. 



Again, the corresponding formula for the wrong insertion is obtained by obvious replace- 
ments. When dimensional regularization is used as IR regulator the on-shell renormal- 
ization constants are equal to unity and the bare matrix elements are zero except for the 
two-loop diagrams with a massive charm quark-loop insertion, 



(1) 



Air 



M 



(2) 
ab 



+ yP 



ab 



+ 0(afi (O b ) 



(o) 



(15) 



This expression assumes that we determined the operator ultraviolet renormalization 
factors Yab from the general expression (JHJ), using another regulator than dimensional 
for the infrared divergences, but then use (fT5l) in the matching relation since the QCD side 
( Hip is calculated with dimensional regularization as infrared regulator. In this relation 
we eliminate the four-flavour coupling by the D-dimensional relation a s = £45 a a , where 
£45 = 1 + 0(a s ) is expanded as in ( fT2l) . Since the SCET operators depend on a variable, 
the renormalization factors Y ao = Y ao (u, u') are functions of two variables and the product 
in (|T5|) is a convolution. 

Just as for four-quark operators, we can write down a matching equation for QCD 
currents to SCET currents. This is convenient since in this way we can organize the cal- 
culation such that the factorizable diagrams are cancelled, at least for the right insertion. 
We consider the relational) 



1 -yC 1 -75)9 = C qq X-2~(l -75)x, 
<Z^+(! - 7s)& = C FF - 75)^ 



(16) 



Since a single collinear sector in SCET is equivalent to full QCD [20] we have C qq = 
1 + 0(a 2 s ) because diagrams with massive bottom quark loops (not contained in SCET) 
arise only at two loops and beyond. The matching coefficient Cff can be determined 
from matching calculations for the b — > u transition [18,25-29]. In the notation of [28] 
we have 



C 



u 



FF 



C\ + — C*2 + C31 



1 - £ % (6L2 + 30L + 7x2 + 72) + ' (17) 



2 The first relation is somewhat schematic, since we suppressed the non-locality of the operator and 
the Wilson line connecting the quark fields on both sides of the equation. 
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with 

L = ln(^). (IS) 



m b 



From ( {TBI) one can construct the following factorized QCD operator 

QCD = [?^(1-7b)?] [g^ + (l-T5)6] = CWQ*Oi, (19) 

which is defined by the product of the two currents given above. By construction the 
matrix element of the factorized QCD operator is the product of a light-cone distribution 
amplitude and the full QCD form factor. In terms of this operator the matrix elements 
of Qi read for the right and wrong insertion, respectively, 



(Qi) = % (0 QCD ) + H ™( a) > (2°) 

0>1 

(Qi) = % (0 QCD ) + H iX (6 x - 1 ) + Y,H ia (O a ) . (21) 

a>l 

Comparing (120]) and fT2~T]) with (J6J and (jUJ) one finally obtains for the hard-scattering 
kernels 

Ti=r^, ^=r%~- ( 22 ) 

^ FF ^qq ^FF ^qq 

3 Computational methods 

The calculation of the two-loop QCD amplitude involves the evaluation of the 62 non- 
factorizable diagrams shown in Figures 15 and 16 of [2]. In addition, there are four 
diagrams which involve a gluon self-energy, and in each of these the self-energy loop 
comprises massless quarks, gluons, ghosts, and the massive charm and bottom quarks. 

We work in dimensional regularization with D = 4 — 2e, where UV and IR (soft and 
collinear) divergences appear as poles of up to the fourth order in e. Since we work in the 
CMM basis we can adopt for 75 consistently the NDR scheme with fully anticommuting 
7s- 

The amplitude of the diagrams is reduced by techniques that have become standard 
in multi-loop calculations. We apply a Passarino - Veltman [30] reduction to the vector 
and tensor integrals. The Dirac and color algebra is then performed by means of an in- 
house Mathematica routine. The dimensionally regularized scalar integrals are further 
reduced to a small set of master integrals using the Laporta algorithm [31,32] based on 
integration-by-parts (IBP) identities [33,34]. To this end we use the package AIR [35]. 

The techniques we apply during the evaluation of the master integrals are manifold. 
The easier integrals can be written in closed form in terms of T-functions and hypergeo- 
metric functions and can subsequently be expanded in e with the package HypExp [36,37]. 
In more complicated cases we derive Mellin-Barnes representations by means of the 
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package AMBRE [38]. We perform the analytic continuation to e = with the pack- 
age MB [39], which is also used for numerical cross-checks. We then apply Barnes' 
lemmas and the theorem of residues to the multiple Mellin-Barnes integrals, and insert 
integral representations of hypergeometric functions as well as -^-functions and Euler's 
5-function where appropriate. As a third technique we apply the method of differential 
equations [40-42] and evaluate the boundary condition with the Mellin-Barnes technique. 
Eventually, the master integrals are evaluated as Laurent series in e. 

The master integrals of all diagrams except the ones in which the gluon self-energy 
contains a charm-quark loop are obtained in a fully analytic form. The coefficient func- 
tions of the Laurent series in e are logarithms, polylogarithms and harmonic polyloga- 
rithms [43] of maximum weight four. Two master integrals of the diagrams with a charm 
loop in the gluon self-energy cannot be displayed in a purely analytic form, see (1751) 
and (1781) . Hence the charm-dependent part of the hard-scattering kernels still contains 
two- fold auxiliary integrals, see (1931) . (1771) . and ([79]) . However, after convolution of the 
hard-scattering functions with the Gegenbauer expansion of the light-cone distribution 
amplitudes, we obtain fully analytic results, including the charm-mass dependence, for 
the topological tree amplitudes a\ and «2- 

The master integrals have been calculated already in [11] and they have been used 
in [12]. We find perfect agreement on the expressions given in [11]. Nevertheless, we 
give in Appendix [A] explicitly the results of a selected subset of master integrals. The 
reason for this is that we found closed forms valid to all orders in e, or that we need the 
structure of the result for the hard-scattering kernels and the topological tree amplitudes 
derived in the subsequent sections. 



4 Master formulas 

From the equations given in Section [2] we can deduce the master formulas for the hard- 
scattering kernels. Expanding them according to 

T^T^ + ^+^yT^ + Oial), (23) 

where a s denotes the five-flavour coupling in the MS scheme, the master formula for the 
hard-scattering kernels for the right insertion reads 

T 

_ 4 (l)nf 7 (1) .(0) 



(0) 

i 


- 4 (0) 


(1) 

t 


4 (l)nf 


(2) 


= A (2)nl 



if = + z m A m + Z f) A <o) + z m A m- + H) Sm m A >m- 

-7-«[ C w+ym-z«]-i:ff«y«. (24) 



6>1 



The actual calculation amounts to evaluating the two-loop on-shell matrix element of the 
transition b(p) — ► qi(p')q2(uq)q 3 (uq) with kinematics p = p' + q, p 2 = mf, p' 2 = q 2 = 0, 



S 



and u (u = 1— it) the momentum fraction of the quark (antiquark) in M 2 . Thus, AjS 
and Hft depend on u, and a term involving Y^ such as Y^ must be interpreted 
as the convolution product J* du' \u') Y^\u', u). 

While the individual terms on the right-hand side of ( 124"1) may be divergent, the 
hard-scattering kernels are free of poles in e. In general, we need tree-level and one-loop 
quantities such as etc. to order O(e) (and even 0(e 2 ) for and T/ 1 ^), 

since they multiply divergent renormalization constants in the expression for the two- 
loop hard-scattering function T { . The are the operator renormalization constants 
from the effective weak Hamiltonian. Their matrix expressions read [44,45] 






(25) 



_£W , 17 ST f n f 26 25 47^/ 

3 ' 9 3 18 6 27 

47^-39 4 ^-f 

ST f n f 79 20T / n / _ 205 5T f n f 1531 2T f n f ± j_ _ J5_ 

9 12 27 18 108 ~ t ~ 288 81 72 384 864 

1( > T f n f 1 83 o T f n f I 119 8 _J5_ _J7_ 

3 4 ° 9 ~ t ~ 16 9 192 72 

where n/ = rty + 1 = 5 is the total number of quark flavours, and Tf = 1/2. The 
row index of the matrices and labels (Qi, Q2, E^\ E% , e[ 2 \ E^), whereas the 
column index stands for (Qi,Q2)- Hence in all terms of the form Z^Aji we also have 
to consider the matrix elements of the evanescent operators, including the factorizable 
diagrams. The terms proportional to Z& \ 8wP-\ and Z^ t account for coupling, mass, 
and wave function renormalization, respectively. The coupling is renormalized in the 
MS scheme with five active quark flavours, whereas the mass and the wave-functions 
are renormalized in the on-shell scheme. The term Cff was already discussed in the 
previous section. The SCET renormalization kernel factorizes according to 

y« ( n , u') = Zf 8{u - u') + Z { bI(u, u) (26) 
into the universal renormalization factor for SCET currents, 











\-L - 5" 


} 




2_ 





(27) 



and the contribution from the ERBL [46,47] kernel, for which we use the convention of 
Eq. (51) in [48]. The last term in (124")) emerges from the fact that the SCET evanescent 
operators 02,3 mrx m to 0\. Renormalizing the evanescent operator matrix elements to 
zero introduces the finite off-diagonal renormalization constants Y^\ As a consequence, 
we also have to calculate the one-loop hard-scattering kernels of the evanescent operators, 
, b > 1, for which only 6 = 2 arises at this order. We find that only the diagram 
with gluon exchange between the n + -collinear quark and antiquark lines contributes to 
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the off-diagonal renormalization constant K;j , as was the case in a related calculation 
performed for the spectator-scattering kernels in [6]. The required convolution results in 

^ (1) (1) 64 L 800 16 7T 2 16 7T 2 32 1n(l-«) 32 In(«) 



b=2,3 

32Li 2 (l-u) 32 Li 2 (u) 32ivr 
9^ 9(1 -u) 9~~ 



(28) 



for z = 1; it vanishes for i — 2. 

The master formula for the expansion of the wrong insertion hard-scattering kernels 
reads 



7 W 7(0) , T(l)f 4 (l)f 7(0) ry(l) v (l)i T(0) 

J J v v ' v v ' 

0(e) 0(e) 



7?(o) 
« 




f(i) 
j 


_ 7(l)nf 


f (2) 


= A (2)nl 



(l)y(l) 
sb 1 61 
b>l 

+(^»+^)[is; ) '-4 i 1 ), i<?] 

- (cS - eff ) US' - a'? 1 - Ki 2) - iff] 4» • (29) 

All terms in the first three lines of T> have a corresponding term in the master formula 



(1241) for the right insertion. The fourth and fifth lines describe Fierz differences of 
factorizable QCD diagrams not contained in the QCD form factor or LCDA. The sixth 
line stems from two-loop SCET matrix elements with a massive charm-loop insertion. 
The terms in the last line of (129]) emerge from the fact that the difference 0\ — 0\ is 
evanescent and hence we must renormalize its matrix element to zero. Similar terms 
already appear in the one-loop expression in which case, however, they are 0(e), 
as indicated, and can be dropped. 

In the two-loop expression, we here find — Ylb>i -^sb = ^> srnce the mixing of 
O2 into 0\ at one-loop is of higher order in e. The terms in the fourth to the last line 
of the expression for lf 2) in (EEL which appear only for the wrong insertion, are given 
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Figure 2: Examples of factorizable diagrams of class A (left), B (center), and C (right). 



explicitly for i = lH as follows: 



T(2)f ,(2)f J(0) 
^11 — /i 21 ^11 



64 




81e 



(^ ) + ^)[^S )f -4f^] 



64 



81 



0(e), 



-[Mg)-M}?]IS? 

— S45 )l r ll r 11 J -^11 




(33) 



0(e), 



(34) 




136 



81 



(35) 



Let us elaborate more on (1301 - (|35|) . In (|30|) we have to calculate the factorizable two- 
loop diagrams both, for the right and the wrong insertion. Since the two are related by 
a Fierz transformation, many terms drop out, resulting in the simple expression ( 130]) . 
There are three classes of two- loop diagrams that contribute to (1301) . They are labelled 
class A, B, and C, respectively, and examples are shown in Figured Diagrams of class A 
do not contribute to the on-shell matrix elements since scaleless integrals vanish in di- 
mensional regularization. For the same reason, only those diagrams of class B contribute 
that contain a massive quark loop insertion in the gluon propagator. Diagrams of class C 
have to be evaluated explicitly, and their right-insertion result can be checked against the 
two- loop heavy-to-light current matching calculation at the kinematic endpoint [26-29] . 
In (|3"T1) the renormalization constants multiply the Fierz difference of the one-loop factor- 
izable diagram, and in (1321) we evaluate the Fierz difference of the one-loop factorizable 
diagram with a mass counterterm insertion on the heavy line; the latter turns out to 
be of 0(e). Eq. ( J33l) comprises the Fierz difference of two-loop on-shell SCET diagrams 
with a charm-loop insertion into a gluon propagator. In order to determine the SCET 
renormalization constants and Y-f± in (l34"l) and (1351) . we must use an IR regulator 
different from the dimensional one. The difference — Y^' must then be a finite, 
regulator-independent constant. However, it turns out that this difference is 0(e), hence 
there is no contribution to from fl34|) . 

3 The corresponding terms for i = 2 are in these cases obtained by removing a factor of Cf = 4/3. 
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The last term (135]) . which does not appear in this form in (129]) . requires further 
explanation. While Yyy is a MS renormalization constant for the physical SCET operator 
Oi, Y-iy is a non- minimal renormalization factor determined by the requirement that the 
matrix element of the evanescent operator 0\ — 0\ vanishes. As a consequence, Yq '— Y] 
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depends on the (non-dimensional) infrared regulator at 0(e 2 ), and so does the finite term 

~ (2) (2) 

of the two- loop renormalization factor — , for example from a diagram such as the 

first one in Figure [2j The infrared regulator dependence of Yq — Y$ enters the 0(e 2 ) 

terms of as can be seen from (1291) . and subsequently the finite term of t!> through 

TpY^. Of course, the final result for must be independent of this regulator. To see 

this we extract the infrared-sensitive 0(e 2 ) term from by decomposing it in the form 
~ ( i) = ~ ( i )F _ _ Y ^ ~ (0 ) CQmbine the subtracted term with _{yW _ yP)] ^4(0) 

from (129|) . This results in the left-hand side of (|35|) . We then calculate the combination 
of renormalization factors that appears in this equation and obtain the infrared-finite 
result given above. For this calculation we used two different methods. One applies the 
method of IR-rearrangement similar to [49,50]. The other combines terms algebraically 
such that no infrared regulator has to be introduced to calculate fl35|) . In this way we 
also verify that the finite term is regulator independent. Summing all terms we find that 

~Y 2) 

Tj is infrared and ultraviolet finite as required by the factorization theorem. 

The final result for the hard-scattering kernels is lengthy and complicated. We there- 
fore relegate the explicit expressions for the right-insertion kernels Tp and T- (2) (z = 1, 2) 
to Appendix IB. ill It turns out that the result of the calculation of the wrong- insert ion 



kernels can be written as linear combination of the right- insert ion ones and a residual 
constant as follows: 

f m = _It (1) - - 
1 3 1 3' 

f 2 (1) = 2Tf, 

i 3 i 9 2 27 ; 27 ; 9 

f 2 (2) = 2 T; (2) + - T 2 (2) + 6 T[ 1] . (36) 
3 

5 Convolution in Gegenbauer moments 

For the topological tree amplitudes to NNLO we adopt the following notation 



i>i 



+ ..., (37) 



4 Note that the appendix contains the one-loop kernels explicitly only through order O(e ), although 
their higher order terms enter the two-loop master formulas. 
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where the ellipsis stands for contributions from hard spectator-scattering which are not 
subject of the present paper. As before a s denotes the five-flavour coupling in the MS 
scheme. The quantities are convolutions of the hard-scattering kernels with the 
light-cone distribution amplitude of a light meson, 



(o) 



V, 



1.7 

(0) 



2.7 



duTf ] <p M {u) , 



du Tf> (f) M (u) 



C 



F y d) 



2N r . 11 



duTf\u) 4> M {u) 



C 



2N r 2j 



du ff\u) 4> M (u) . 



(38) 



The light-cone distribution amplitude of a light meson, 0m, is expanded into the eigen- 
functions of the one-loop renormalization kernel, 



4>m(u) = 6u(l - u) 



n=l 



(39) 



where a^f and Cn^ 2 ^ (x) are the Gegenbauer moments and polynomials, respectively, and 



the integrals (I3"8l are then performed term by term. Below we provide integrated results 
truncating the Gegenbauer expansion (|39~I) at n = 2, which is sufficient in practice. The 
light-cone distribution amplitude and Gegenbauer moments are scale-dependent, and our 
notation is such that refers to the scale /i that is also the scale of a s . Since 0m(w) is 
normalized to unity and the tree-level scattering kernels are constant, we simply have 



(0) 



ii 



T/(0) 
v 12 



0. 

1. 



V, 



(0) 



21 



y(0) 



4 

9 ' 
1 

3 ' 



(40) 



At the one-loop level, one obtains 



(i) 



n 



y0) 



12 



45 

6L — 3in 

2 



H o- \ M 21 M 

Stir \ a-, an i , 

2 1 20 2 ' 



V, 



(i) 



21 



(l) 



11 



V, 



(1) 



22 



2V«. 



(41) 
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Note that from one-loop onwards, the constants in the Vij are scheme dependent and 
hence look different in the CMM basis compared to those in the traditional basis [16] em- 
ployed in the earlier QCD factorization calculations. This explains the different constant 
in v£i compared to the corresponding expression V M in [3]. 

At two loops the expressions get more complicated. In order not to spoil the simple 



(2) 

structure of the result we split the according to 

+ T f AV, 



(2) _ T/ (2), , m »,,(!), 



v> ! = V 

V ~ I'J 



(42) 



The first term on the right-hand side describes the result for a massless charm quark. 
The second term is the difference between the contribution of a quark of mass m c and a 
massless quark, see also (181]) . This construction ensures that we set ni = 4 irrespective 
of whether we treat the charm quark as massive or massless. Explicitly, we have 



v, 



(2), 
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-39L 2 - 234L + 8194C(5) - 2028vr 2 C(3) 

172597 



100862C(3) 



15 



+ 44tt 2 ln(2) 



55vr 4 43909tt 2 
+ 



216 



+T f 4L 2 + 



70L 



+ 



180 
56vr 4 



+ niT f [ AL A + 



70L 4tt 2 493 



3 + 9 



+ m 



-39L 



3 75 3 9 

f-VS-f^MJ + tf)) 

20C(3) 13017T 2 



164vr 2 1321 4128 >/oN 
+ — + — =— C(3) 



270 



1591 m / 44 



+T f ( AL + - 184^5 ln(f + f) - 288 ln 2 (i + 



yl5\ 

2 > 



192 



C(3) 



1639L 



50964C(5) + 12618tt 2 C(3) + 



208846C(3) 



- 324vr 2 ln(2) 



1163vr 4 5411297T 2 3604333 
+ + — + n x T f 



120 



22L 
IT 

716 



+ 



36 
19532 



1080 



22L 



+ 316tt' 



544tt 4 
25 



3 

29856 



-4tt z 



80 

T 



C(3) 



+ Z7T 



149L 



+ 



5 

1037T 2 



7T 2 \/5 



10 



5923 



^7r 2 ln(i + f)^) 

+ ri z 7) (4L + 14) + 7) ^4L + 990 



576 



C(3) 



-1432^5 ln(| + % + 2976 ln 2 (± + 



14 



i 1169L w . 2>/ , 4896383C(3) 2392 2 . . 
f af <| — + 170274C(5) - 42138tt 2 C(3) + — tt 2 ln(2) 

199tt 4 13604203tt 2 733820723 m /7L 2 8059 \ 

+ + nj T f 2vr 2 + 

12 270 75600 ' V 5 300 / 

_ /7L 1585223 3322 2 1112 4 36288^ , 



+ i"7T 



5 

953tt 2 47929 6 m _ / 39182 1152 w , 



210 720 



+5536^111(1 + ^) - 13248 ln 2 (± + ^ | . 



(43) 



The corresponding charm-mass dependent part reads 

A\/A 2) | = ^-^z 2 q 3 (z) + 16z 2 \n{z)q 4 {z) + 24z 2 q 5 (z)-^{22z+l)q 6 (z) 

+l8z(Az - l)q 7 (z) - 4v^ (28,2 - l)q 8 (z) + 16^ 2 g 10 (^) + 2z 2 q u (z) 
-16z 3/2 q 12 (z) + 92z - 2n 2 ^23z 2 + 3z + y/z-^j- 16v^Li 2 (l - yfz) 
+2 (3z 2 -3z-l) ln 2 (z) + 76z \n(z) + 4 (3z 2 -3z + y/z + l) Li 2 (l - z) 
+ m 2Az 2 (2 \n(z) - 3) ln 2 (??) - 32z 2 gi(^) + 296,2 

-4\n(z) + 4(22z + l)q 2 (z) 

+ of | - 288,2 2 ln 2 (?7)Li 2 (?7) - 4^ 2 g 3 (z) + 8,2 2 (l2z 2 + 16,2 + 3 + 6 ln(z)) q 4 (; 

-- (72z 3 + 8Az 2 + 22z + l) q 6 (z) - — (45z 3 - 5,2 2 + 3) Li 2 (l - z) 

O LOZ 

2 ( 467 1\ 

+- 36z 4 - 312,2 3 + z 2 -Z + -) q 7 (z) + 576,2 3 - 336^ 2 

z \ 3 5/ 

4 / 4 , 2 3\ , . 18892 

216,2 4 - 504,2 3 - 44z 2 -32z+-) q 8 (z) + 



3 V 'T V """ 5 J ' ' 3 

+8z (3z 3 - 8z 2 -6z-3) q w (z) - I2z 2 q n (z) + 112,2 3/2 gi 2 (,2) 



15 



15^ 



(20/ + 20z + 3) (8Li 2 (l - y/zj - 2Li 2 (l - z) + tt 2 ) 



2tt 2 

+ (2325/ + 5/ + 30z - 3) + 4z (144/ - 72z + 613) \n(z) 

15z 

(45/ -5/ + 3)ln 2 (^) 



15z 



+ m 



288/ + 24/ (12/ + I62 + 6 ln(» + 3) ln 2 (r?) - 96zV(*) + 360/ 



+328^ - 4 ln(z) + 4 (72/ + 84/ + 22z + l) g 2 (z) 



+ af j- 8/ (-31n 2 (r/) - 12Li 2 (ry) + 2tt 2 ) ; 



-16/ (40/ + 30/ - 6 ln(s) - l) q A (z) - -tt 2 z(1Q49z + 9) 



+ m 



+ -z (120/ +70z 2 -llz - 6) g 6 (s) +6(z - l)z (ln 2 (z) + 2Li 2 (l - z)) 

-2z (240/ - 330/ - 971z + 9) q 7 (z) 

+ 16/ /2 (120/ - 185/ - 51z - 147) q 8 (z) 

-4/ (40/ - 75/ + 11) q w (z) + A2z 2 q u (z) - 9Qz z/2 q l2 (z) 

—z (2880/ - 1140/ - 823z - 54) ln(z) 
3 

~z (17280/ - 8280/ - 2861z + 1908) 
9 

- 48/ (40/ + 30/ - 6 InO) - l) ln 2 (r?) - 192z 2 qx(z) 



--z (720/ + 480/ - 37,2 - 183) 
o 

-I62 (120/ + 70/ - llz - 6) q 2 {z 



(44) 



with z = m 2 /m 2 . The variable r\ and the functions qi(z) are defined in Appendix IB. 21 
One can check that AV^ \ vanishes for z — > 0, and that for z — > 1 it approaches the 

corresponding contribution of a bottom quark contained in (|43l) . 

(2) 

The case of V 12 is simpler since the corresponding one-loop expression vanishes due 
to the colour structure, and therefore one also has 



AV, 



(2). 



12 



0. 



(45) 
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Hence there is no need to distinguish between the massive and massless charm contribu- 
tion. We obtain 



(2) 



1:2 



18L 2 + 156L + 3744C(5) - 936tt 2 C(3) - 14833 ^ 3 ) + 72tt 2 ln(2) + 

5 90 



+ 



12487vr 2 5347 



12 



60 



+ iir 



, . 47vr 2 1333 
18L - 16C(3) + + 



45 



12 



+ a x 



AI 



-33L - 19008C(5) + 4752tt 2 C(3) + 77157 ^ _ 2 4vr 2 ln(2) - ^ 

O -LU 



21807tt 2 4568 

+ — rzr- + m 



4 



15 



r 23vr 2 15 

18L + + — 

5 4 



63L 



+ 74304C(5) - 18576tt 2 C(3) - 



2236872C(3) 



10 

204218tt 2 



32369221 
12600 



+ m 



187T 2 

^5" 



35 

173 
~30 



2064vr 2 ln(2) + 



797tt 4 
10 

(46) 



The integrated two-loop kernels for the colour-suppressed tree amplitude follow from 
(I3"6"j) . and read 



v 21 



v 22 



1 



,(2) 



10 



10 



49 

y 



,(2) 



(1) 



(47) 



Our results agree with those in [12,13]. The comparison we performed is of analytic 
nature, except for those contributions to the real part of the amplitudes that stem from 
charm and bottom loop insertions into a gluon propagator, for which the result has been 
obtained in [13] only in numerical form. For these terms the numerical agreement for the 
values of the quark masses covered in Table 1 of [13] with our analytic result is within 
three permille. 



6 Phenomenological applications 
6.1 Input parameters 

In this section we begin the discussion of the numerical evaluation of the topological tree 
amplitudes a±, ct2, and the tree-dominated charmless B decays into final states 7r7r, 7rp, 
and pp. The main result of this work, the NNLO vertex correction, depends on very few 
input parameters: the strong coupling a s via A™.; the bottom pole mass m& or MS mass 
m&(m&); the renormalization scale p = pb (to distinguish it from the hard-collinear scale 
that appears in spectator scattering); the Gegenbauer moments of the light-meson light- 
cone distribution amplitude aff , the only hadronic parameters; and the electroweak-scale 
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Parameter 


Value /Range 


Parameter 


Value /Range 


A® 

MS 


0.225 


/^hc 


1.5 ±0.6 


m c 


1.3 ±0.2 




0.195 ±0.015 


m s (2 GeV) 


0.09 ±0.02 


U 


0.131 


(m u + m d )/m s 


0.0826 


/*(0) 


0.25 ±0.05+ 


m b 


4.8 




0.209 


rh b {m b ) 


4.2 


< p (o) 


0.30±0.05 t 


\v*\ 


0.0415 ±0.0010 


A B (1 GeV) 


0.35 ± 0.15+ 


1 T 7" / T T 1 

\v ub /v cb \ 


0.09 ± 0.02 


(Ti(l GeV) 


1.5 ± 1 


7 


(70 ±10)° 


(7 2 (1 GeV) 


3 ±2 


r(B") 


1.64 ps 


4(2 GeV) 


0.2 ±0.15 


r(B d ) 


1.53 ps 


a p 2 (2 GeV) 


0.1 ±0.15 




4 8+ 4 - 8 


a£ ± (2 GeV) 


0.1 ±0.15 



t Value and range changed to /f ^(0) = 0.23 ± 0.03, A^ p (0) = 0.28 ± 0.03, 
A B (lGeV) = (0.20±g;gg) GeV in parameter set "Theory II". See text in 
Section [631 



Table 1: List of input parameters. Dimensionful parameters are given in units of 1 GeV. 



initial conditions for the Wilson coefficients of the effective Hamiltonian. Their values 
and uncertainties are specified in Table dj except for the electroweak-scale parameters 
M w = 80.4 GeV, m t {m t ) = 167 GeV and sin 2 ^ = 0.23. When we evaluate the full 
topological tree amplitudes and branching fractions including penguin amplitudes (using 
results from [9]), we need further Standard Model and hadronic parameters, which are 
also listed in the Table. The values and uncertainties follow [9, 51] with some minor 
adjustments that account for recent reevaluations or progress in error estimates of 7, 
f B [52,53], (0) [54], and a£ (JL) [55]. 

Together with the NNLO non-leptonic matrix elements, we now use three-loop evo- 
lution of the strong coupling, as well as the next-to-next-to-leading-logarithmic (NNLL) 
approximation to the Wilson coefficients appearing in the effective weak interaction 
Hamiltonian in the CMM basis. This approximation is constructed from the two-loop 
initial conditions [45,56], three- loop anomalous dimension matrix [44,45] and the NNLL 
solution to the renormalization group equation [57]. The penguin amplitudes, which 
are not yet known at NNLO, are evaluated exactly as in [9], whenever needed, except 
that we employ the same three-loop evolved strong coupling in the entire numerical pro- 
gram. The observables we discuss in this paper are chosen such that they do not depend 
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significantly on the penguin amplitudes. 



6.2 Topological tree amplitudes at NNLO 
6.2.1 Two-loop vertex correction 

We first consider the loop expansion of the pieces of the topological tree amplitudes not 
related to spectator scattering, as defined in (j3"7|) . To acquire an idea of the convergence 
of the perturbative expansion we provide a numerical representation of 



/ duT j {u)<p M {u) = vif + 



4tt/ 2N r lj 



(48) 



up to 0{a 2 g), and similarly for Tj(u) and ■ Choosing /i = m& we get 



Jo 



duT^u) (j) M {u) = + [-5 - 2.094 i + (1.222 - 2.094 i)af - 0.233 af] ( 

+ [-38.510 - 75.389 i + (142.014 - 24.416 z) af 
+ (-18.588 - 13.528 i) af + A mc ] (J^j ' , 

duT 2 {u) (f) M (u) = 1+ [71.598 + 71.320 i + (-47.327 + 34.313 z) af 

+ (0.848- 7.569 z')afl f— 

V47T 

dufi(u)^jif(u) = ^ + [0.333 + 0.698 i + (-0.407 + 0.698 z) af + 0.078 a 

+ 37.362 + 56.828 z + (-68.372 + 23.389 i) af 



V47T 



](£) 



+ (6.573 + 1.145 z) af - - A 



47T 



/' 

./o 



duT 2 (u) 4> M {u) 



~+ [-10- 4.189 z+ (2.444-4.189 
3 



z a. 



0.467 a 



2 J Utt/ 



+ [-83.154 - 139.571 z + (275.585 - 49.960 z) c 
+ (-38.293 - 29.580 z) af + 2 Am c ] ^) 2 . 



(49) 



We observe that the NNLO terms have relatively large coefficients, but since a s (mb) w 
0.22, the perturbative expansion is still under control for most of the terms. The coef- 
ficients of the second Gegenbauer moment are generally rather small, contrary to those 
of the first one, as we will discuss in more detail at the end of this subsection. 
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The numbers in ( |4"9]) apply when we treat the charm quark as massless. The charm 
mass effect, the difference between the contribution of a quark of mass m c and that of a 
massless quark, is parameterized by A mc , for which we find the fitting formula 

A mc = -0.9661 - 2.0849 i- (11.4437 +12.9458 i)(z- 0.075) 

+ (28.0366 + 64.7300 i)(z - 0.075) 2 
+ [4.5919 - 1.9379* + (14.4270 - 12.2979i)(z - 0.075) 

- (96.8420 - 59.5511 i)(z - 0.075) 2 ] of 
+ [-0.3599 - 0.3833 i - (2.2574 + 0.5738 i) (z - 0.075) 

+ (14.1412 + 7.0945i)(z - 0.075) 2 ] af (50) 

in the region of z = m 2 c /m\ of interest. The second-order polynomials in z were obtained 
from a least-squares fit to the exact functions, which works better than 0.4% for z G 
[0.05,0.10]. The size of the charm-mass dependent terms is moderate. Most of the 
numbers in (1491) receive small corrections, but the shifts induced by the z-dependent 
terms can be up to 11% for single entries. 

The successive LO, NLO and NNLO approximation^] to the topological tree ampli- 
tude parameters ai(7T7r) and a 2 (7T7r) (vertex correction only, i.e. no spectator scattering) 
are shown in Figure [33 depending on the renormalization scale \x. Relative to the NLO 
correction the NNLO term turns out to be sizeable. For the colour- allowed amplitude 
a;i(7r7r) this can be explained by the observation that the NLO correction is colour- 
suppressed, while the NNLO correction is no longer colour-suppressed for any of the two 
amplitudes. As in any perturbative QCD calculation the higher-order corrections are 
expected to reduce the scale uncertainties compared to lower-order approximations. We 
observe that this is indeed the case after taking into account the NNLO corrections for 
the real parts of the amplitudes, but the reduction is absent for the imaginary parts. We 
attribute this to the fact that the two-loop contribution to the imaginary part is really 
a next-to-leading order effect since the imaginary part is zero at tree-level, and to the 
large size of the two-loop correction. 

Besides the renormalization scale the vertex correction to the topological tree ampli- 
tudes depends only on the Gegenbauer moments. At leading order there is no dependence 
on these moments, whereas at NLO this dependence is given by 

[ai] v = 1.040 + 0.013z - (0.007 - 0.013i) af 2 + 0.001 af 2 , 

[a 2 ] v = 0.029 - 0.079z + (0.046 - 0.079z) af 2 - 0.009 af 2 , (51) 

5 The LO and NLO result is computed in the traditional operator basis for the effective weak Hamil- 
tonian, so that it corresponds to the expressions used in [1,3], but is slightly different from the NLO 
truncations of the NNLO results of the present calculation. 
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Figure 3: Dependence of the topological tree amplitudes 0:1(7771) and 02(7177) on the hard 
scale n (vertex corrections only). The dotted, dashed and solid lines refer to the theoretical 
predictions at LO, NLO and NNLO, respectively. 

where we choose again \x = nit,. At NNLO we obtain 

[ai] v = 1.057 + 0.038i - (0.032 - 0.022i) of 2 + (0.003 - O.OOlz) af 2 , 

[a 2 ) v = 0.013 - 0.126z + (0.139 - 0.096?;) af 2 - (0.021 + 0.009?;) af 2 . (52) 

Since the Gegenbauer moments are typically smaller than 0.3, we conclude that the 
dependence on the second moment af 2 is very small for both tree amplitudes. This 
implies that there is practically no dependence on the shape of the pion or rho meson 
light-cone distribution amplitude in the vertex term of the factorization formula (pQ). 
This is no longer true once spectator scattering is included, especially for the colour- 
suppressed tree amplitude a 2 , as discussed in [6]. The first Gegenbauer moment is more 
important and has a significant effect on the value of the colour-suppressed amplitude 
[aiz] v, which is furthermore enhanced at NNLO. This is a source of non-factorizable 
SU(3) flavour symmetry breaking (not related to decay constants and heavy-to-light 
form factors) in the relation of amplitudes with pions and kaons. 
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6.2.2 Complete NNLO amplitude including spectator scattering 



We now proceed with our investigation of the topological tree amplitudes by adding the 
corrections from spectator scattering on top of those from the vertex corrections. Our 
numerical results for the irn final states read 

a^vrvr) = 1.009+ [0.023+ 0.010 z] NLO + [0.026 + 0.028 i] NNLO 
r, 



sp 



0.445 



[0-014] LOsp + [0.034 + 0.027* 



NLOsp 



+ [coos; 



tw3 



l.ooo; 



-0.029 
-0.069 



(0.0111°; 



0.023\ 
050 J 



(53) 



a 2 (7i7i) = 0.220 - [0.179 + 0.077 % 



NLO 



[0.031 + 0.050 i] 



NNLO 



sp 



[0-114] LO sp + [0-049 + 0.051* ] NLOsp + [0.067] 



.0.445 

0.2401S + (-0.077^ 15 



tw3 



-0.078 



(54) 



Here the "NNLO" term in the first line of each of the two equations corresponds to 
the numerical evaluation of the calculation of this paper, while the second line accounts 
for spectator scattering at 0(a s ) ("LOsp") and 0(a 2 s ) ("NLOsp") and a certain power 
correction ("tw3"). The theoretical calculation for spectator scattering is taken from [6], 
but the numerical values differ, since we now use a different operator basis, which redis- 
tributes the perturbative expansion, and slightly different input parameters. The third 
line sums all contributions and provides an estimate of the theoretical uncertainties from 
the parameter variations detailed in Table [TJ By comparing the theoretical uncertainty 
in the full expressions above to the one from the vertex correction alone, see Figure [31 
we see that it arises primarily from spectator scattering. The main contributors to the 
uncertainty are the parameter combination 



9UL 



B 



sp 



mJr(0)A 



B 



(55) 



which appears as an overall normalization factor of the spectator-scattering term, the 
second Gegenbauer moment a2(2GeV), and "tw3" power correction. 

Relative to the full amplitude the two-loop vertex correction is only a few percent 
for the colour- allowed amplitude an, but it is quite significant for a 2 . The negative 
(10 — 15)% correction to the real part decreases the branching fractions of the colour- 
suppressed decays by about 25% relative to previous results that included the 0(ot 2 s ) 
correction only in spectator scattering. The magnitude of the correction to the imaginary 
part of c*2 at two loops still reaches 25% of the leading-order real part. However, we 
also observe that the two 0(a 2 ) corrections, "NNLO" in the vertex term and "NLOsp" 
in spectator scattering, cancel to a large extent in both, their real and imaginary parts. 
This is somewhat unfortunate, since an enhancement rather than a cancellation in the 
colour-suppressed tree amplitude might have been welcome in view of the trend indicated 
by experimental data, as will be seen below. 
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The colour- suppressed tree amplitude 02(^71") exhibits an interesting structure. It 
starts out with a positive real value 0.220 that corresponds to the naive factorization 
approximation. After adding the perturbative corrections to the four-quark vertex (the 
first line of (|54|) ). it is found to be almost purely imaginary, 0.01 — 0.13z. However, 
the spectator-scattering mechanism regenerates a real part of roughly the original size 
and cancels part of the strong phase. The net result of this is that the colour-suppressed 
tree amplitude can become sizeable in QCD factorization when r sp is large, but since this 
enhances the cancellation of the imaginary part, one cannot have both, a large magnitude 
and a large strong phase. In comparison, the colour-allowed tree amplitude ai(7T7r) is 
rather stable against radiative corrections, and never deviates by a large amount from 
its naive-factorization estimate. 

6.3 Parameter dependence of branching fraction ratios 
6.3.1 Factorization test 

In this subsection we consider a non-leptonic to semi-leptonic decay ratio that provides 
direct access to the magnitude of the topological tree amplitudes. This analysis has first 
been performed in [5] at NLO and can now be repeated at NNLO, and with improved 
experimental data (see also [15]). 

Within the approximation that all electroweak contributions are neglected since they 
are never CKM-enhanced and formally of order 0(a em ), the process B~ — ► 7t~tt° is a 
pure tree decay. The corresponding decay amplitude in the SM can be written as [5] 



(56) 



with A ww = i rngf+^^fn, and Ap = V p ijV* d . A substantial uncertainty in the predic- 
tion of the branching ratio arises from the overall normalization due to \V u i,\ and f?*(0). 
These two quantities cancel out in the ratio [5,58] 



7T 7T 



Rn = -=77*T ,1 = 37r ft \ V udY Mtttt) + a 2 (7r7r)| 2 , (57) 



dT(B° -> n+l-u)/dq 



q 2 =0 



and a direct measurement of the tree amplitude coefficients can be performed if the 
differential semi-leptonic B — > irl v rate is measured near q 2 = 0. The ratio in (IBTj) 
therefore constitutes a test of QCD factorization, and furthermore underlines the im- 
portance of measuring the semi-leptonic decay spectrum for understanding the pattern 
of non-leptonic B decays. Replacing pions by longitudinally polarized rho mesons in 
( 1ST)) , we define an analogous ratio R p , which allows us to test factorization in B — > VV 
decays. 

This test can currently be carried out for pions. The latest experimental data on the 
B~ — ► 7r~7T° branching fraction reads [59-62] 

Br (B- -> tt-tt ) = (5.59+°;^) ■ lO" 6 . (58) 
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Since the theoretical expression for the differential semi-leptonic spectrum at the end- 
point assumes the form 

(59) 

the measurement of the left-hand side is equivalent to the experimental determination 
°f \Vub \ 7+^(0) from the full lepton invariant mass spectrum in B —>■ ixl v decay together 
with an extrapolation to q 2 = based on a form-factor parameterization. Several ex- 
tractions of \V u b \ /^""(O) [63-65] employing different parameterizations of the transition 
form factor [66-69] can be found in the literature, which all yield similar results. We 
shall adopt the value 

\V ub \f^(0) = (9.1±0.7)-10~ 4 . (60) 

Combining this and the experimental data on Br(B~ — ► 7r~7r°) given by fl58|) allows us 
to extract the following experimental constraint on R^, 

K]exp= 0.81 ±0.14, 

which can be translated into a constraint on the tree amplitudes, 

[ |ai(7T7r) + « 2 (7T7r)| ] cxp = 1.29 ± 0.11 . (62) 

On the other hand, varying the parameters in the range specified in Table [1] and adding 
the errors in quadrature, our theoretical prediction at NNLO is given by 

| ai (7nr)+a 2 (7r7r)| = 1.24lg;l«, (63) 

which is in good agreement with the experimental data and hence provides support for 
the factorization assumption. The theoretical uncertainty is still large, but much of it 
arises from Xb, Jb arid f^(0), on which the theoretical result depends only through 
the overall factor r sp in (1551) . In Figure @] we therefore plot the dependence of |ai(7T7r) + 
<^2 (vr7r) I on the first inverse moment of the S-meson distribution amplitude, Xb, and the 
remaining theoretical error (grey band). For comparison we show the LO (flat, dotted 
line), NLO (dashed) and NNLO calculation, which makes it clear that the combined 
NNLO correction is rather small. The overlaid horizontal band refers to the experimental 
value (I6"2"j) . displaying again the good agreement with the theoretical calculation. The 
errors are too large to turn this into a determination of the parameter Xb, which may then 
serve as an input for the remaining non-leptonic observables. The central experimental 
value (thick horizontal line) alone would allow the range Xb G [150,400] MeV, which is 
on the lower side of theoretical expectations [1,3,70-72]. 

In the future it should be of interest to extend the factorization test to other final 
states. In Table [2] we summarize our predictions for pions and longitudinal rho mesons. 
Here "Theory I" corresponds to the input parameters of Table [TJ As will become clear 
below, the analysis of non-leptonic decay data suggests the hypothesis of smaller form 
factors and Xb- This leads us to choose a second parameter set, to be discussed in 
Section IS~4l and labelled "Theory II". We consider the corresponding entries in Table [2] 
to be our reference predictions. 
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dT(B i 
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l-u)/dq 2 



m 



B 



192tt 2 



|K 6 |/r(0) 



(61) 
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Figure 4: Dependence of the theoretical prediction of |ai(7T7r) + a 2 ( 7r7r )| (equivalent to 
R n ) on \ B and comparison to [ |ai(7T7r) + a 2 (mr)\ ] cxp = 1.29 ± 0.11. The dotted, dashed 
and solid lines denote the LO, NLO and NNLO predictions, respectively. The dot-dashed 
lines give the error band upon varying the input parameters within the ranges specified 
in Table [H excluding the parameters contained in r sp . The experimental data within la 
error is also shown (horizontal band). 



6.3.2 Ratios involving the colour-suppressed tree amplitude 

Next we consider a number of branching fraction ratios to further highlight the impor- 
tance of hard spectator-scattering, and thereby the dependence on Ab and r sp . Some of 
these ratios have been considered before, for instance in [5,15,73], for a different purpose. 
We define 



= Z — ; — ; r , -KqO = ^ 



opp _ 2 -+ Plpl) PP _ T(B° -> plpl) 



2T(B° ^7T°p°) 





Theory I 


Theory II 


Experiment 


R P 


u -'°-o.ii 

i 7 c:+0.37 
x - ' d -0.24 


n Q4+ - 23 


0.81 ±0.14 



Table 2: Theoretical results for the factorization test ratios fl57|) for pions and longitudinal 
rho mesons. "Theory II" refers to our preferred input parameter set, see text. 
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Here, and throughout the paper, T denotes the CP-averaged decay rate, i.e. T(B — > f) 
is in fact the average of T(B — > f) and the decay rate T(B — >• /) of the CP-conjugate 
mode. The Rqq ratios are approximately proportional to the square of the ratio of the 
colour-suppressed to the colour-allowed tree amplitude. This approximation is very good 
for the modes involving the rho meson final states, where the penguin contribution is 
known to be small J^| The ratios involve \ai + a^l, similar to R n and R p in the 
previous subsection, but divided by |cux| and squared. These ratios also display some 
sensitivity to the CKM angle 7 for tttt, but we do not discuss this here, since it arises 
through interference with the penguin amplitudes. 

Our results for these ratios are displayed in Figure Again the grey band refers 
to the theoretical prediction (with uncertainty) for the value of on the horizontal 
axis, while the horizontal band gives the current measurement [59-62,74-88]. It is clear 
that naive factorization (the dotted line) fails to describe the data. In general, QCD 
factorization provides a good description of data, especially if A# is around 200 MeV, 
where, however, the theoretical uncertainty becomes large. In the QCD factorization 
approach small A b and a large colour-suppressed amplitude is connected with an impor- 
tant role of spectator scattering in the dynamics underlying this amplitude as already 
noted in [5] . An exception to the good agreement between theory and data is -Rqo which 
can be accommodated even for Xb — 200 MeV only marginally within errors. Thus, it 
is plausible or even likely that power-suppressed corrections are also important for the 
colour-suppressed amplitude, especially if it is accompanied by a large phase, an issue 
that we do not discuss in the present paper, since it requires the consideration of pen- 
guin amplitudes. Explaining why this effect should be large for pions but small for rho 
mesons (rather than fitting it with new parameters [89,90]) appears to be a theoretical 
challenge. 



6.3.3 Ratios involving other np final states 

Next we discuss two ratios that involve only the charged 7r ± p =F final states. These final 
states are dominated by the two colour-allowed tree amplitudes 01(7179) and ai(pir). The 
ratios we consider have few theoretical uncertainties other than the one from the heavy- 
to-light form factors [5]. Specifically, we define 

fla s m° =i W> • AC s 5 [c{w ' p+) " 0(71 V)1 • (65) 

where C(f) denotes the coefficient of the cos(Amet) term in the time-dependent decay 
amplitude into the final state /, see [5] for the relevant definitions. The dependence of 
the theoretical calculation on the form factor ratio F^p = A^ p (0)//+ ir (0) is shown in 
Figure E] together with the present experimental data [59,86-88,91,92]. The agreement 



6 For the vector-vector final state, we consider the longitudinally polarized final state. Accordingly, 
the experimental unpolarizcd branching fraction is multiplied by the measured longitudinal polarization 
fraction. 
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Figure 5: Constraints on Ab from the ratios ( |64|) of B — > nn, B — > irp, and B — > pLpL 
decays. Line styles have the same meaning as in Figure HI 
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Figure 6: Dependence of R 3 and AC on the ratio A^ p (0)/ f^{0) of B -meson form factors 
for fixed f+ n (0) = 0.25. Our default input value for the form factor ratio is 1.20. Lines 
and bands have the same meaning as in Figure HI 

is excellent for F^ p 6 [1.0, 1.2], which is consistent with the input value 1.2 that follows 
from Table [U 

6.4 Final results 

In this subsection we give our final theoretical predictions for the CP-averaged branching 
fractions and a number of ratios of CP-averaged 7r7r, up and PlPl decay rates. In addition 
to the ones already discussed in the previous subsections, we further consider the ratios 

K l = 777777^ 7 77 > K 2 = 



T(B° - 


> 7T + p ) 


T{B° - 


> 7T + 7T - ) ' 


2r(5- 




T(B° - 


-> 7T-p+) 


2T(B- 




T{B° - 


-> TT + p") 



2r( J B° -> 7T + 7T-) 

^4 = V/fn " 1 > ( 66 ) 

- _ r(B^,y) + r(f^,y) 

R3 has already been defined in (165j) . Moreover, we define the following ratios of the 
colour-suppressed decays, 

_ r(g° -> ttV) _ r(g° -> P °p°) 

^ " r( so _ ^o) > ^ - r(j g _> ^o) > t b 

to eliminate the input parameter correlations among the three decay modes. 

The theoretical results and experimental measurements of the CP-averaged branching 
fractions of the twelve final states composed of pions and rho mesons, as well their ratios 
as defined above are summarized in Table [3j The column labelled "Theory I" , on which 
we focus first, uses the input values and uncertainties defined in Table [TJ The first 
error comprises the uncertainties from CKM parameters (7, \V c b\, \V u b\ - see below), the 
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Experiment 
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.08 


B°- 
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+0.06+1.45 
-0.06-0.84 

+0.86+1.22 
-0.69-0.97 

+0.11+0.42 
-0.08-0.17 

+0.42+2.71 
-0.41-1.56 



+0.90- 
-0.77- 



■2.18 
1.41 



-0.56-2.10 

+0.39+1.37 
-0.39-1.42 

+0.27+3.82 
-0.19-3.50 

+0.04+1.11 
-0.03-0.43 



(*) 
(*) 

(**) 

(*) 

(*) 
(**) 

(t) 



cr 09 +0.07+1.42 
d - oz -0.06 -1.35 

r 7n +0.70 +1.16 
°- ' u -0.55 -0.97 

n fi o +0.12+0.64 
u - uo -0.10 -0.42 

O C4 +0.41 +2.54 
-0.40 -2.52 

1 9 1 q +0.85 +2.23 
1Z1 ° -0.73-2.17 

1 q 7fi +0.49 +1.77 
10 -' u -0.44 -2.18 

1^+0.34+1.35 

0- 1 ^ -0.33 -1.49 

91 nn+ - 20 + 3 - 06 
Zi ' 3U -0.12 -3.55 

1 4 q +0.07+1.77 
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(**) 
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1 nr+0.05+1.62 
i - uo -0.04-1.04 



22.8±iJ 
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1 q 9 +0.02+0.44 
i - oz -0.03-0.27 
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U.UO„ .00-0.03 



1.91 
0.22 



+0.18+0.72 
-0.20-0.64 

+0.06+0.28 
0.05-0.16 



1 79+0.03+0.64 
x - ' z -0.03-0.53 

n 1 n +o.oi+o.i9 

u - ±u -0.01-0.11 



2.02 ±0.17 
0.60 ±0.08 

1.8018:1 

0.05 ±0.02 
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+0.00+0.09 
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0.17 ±0.05 
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2.411°; 
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16+0.32 

18- 0.37 
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19- 0.41 

13+1.12 
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2.411° 
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0.64 



0.531°; 
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19+0.42 
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13+0.72 
12-0.59 
00+0.84 
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05-0.41 
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3.04 ±0.37 
2.23 ±0.24 
2.15 ±0.43 
1.12" 



h0.46 
'-0.48 



U.OU_ 23 

0.49 ± 0.08 



R p c 



n fi/ i+0.22+0.64 
u -°^-0.17-0.37 

n 7/1+0.10+0.58 
u - '^-0.09-0.46 



n a O+0.09+0.28 
u -^ z -0.08-0.16 

n 7n +0.06+0.46 
u - ' u -0. 06-0.39 



0.78 ±0.22 
27+°- 13 
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Table 3: CP-averaged branching fractions in units of 10~ 6 [59-62, 74-88, 93-96] and 
various ratios of tree-dominated B — > nn, np and PlPl decays. The first error on a 
quantity comes from the CKM parameters, while the second one stems from all other 
parameters added in quadrature [51]. We consider "Theory II" as our reference values, 
see text. 
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second combines all other uncertainties (scale, hadronic parameters, power correction 
parameters) in quadrature. For the first part of the table containing the branching 
fractions, we apply the following modified procedure. Since the tree-dominated modes are 
dominated by amplitudes proportional to V u b, there is a large normalization uncertainty 
from iKtbl 2 , which we do not include in the table. Instead, we calculate the theoretical 
uncertainty of the quantity BrAv(£? — > f)/\V u b\ 2 ^ The error from \V u b\ can easily be 
restored by assuming that BrAv(-B -> /) oc \ V ub \ 2 . In the same way, one can also rescale 
the branching fractions to account for a value of |V^&/Vd>| different from 0.09. We handle 
the dependence on the heavy-to-light form factors in a similar way. In this case one 
notes that the colour-suppressed final states 7r°7r°, 7r°p°, p°p° are almost independent of 
the form factors, since they are dominated by spectator scattering. The other modes, 
however, are nearly proportional to the square of the form factor. To remove the trivial 
dependence on the form factor uncertainty, which might be reducible in the future, 
we calculate the theoretical uncertainty of the quantity BrAv(i? — > /)//+ 7r (0) 2 for the 
modes marked in the table with (*), of BrAv(5 — > /) / 'A^ '(0) 2 for those marked with 
(**), and of BrAv(5 -> f)/(f+*(0)A$ p (0)) for the mode B -> vr ± p T marked with (f). 
Once again, the full form-factor error can be restored by assuming the dependence on 
the form factor as divided out above, and the branching fractions can be approximately 
rescaled to other form-factor values by multiplying the appropriate factor. 

Scanning the numbers in the "Theory I" column of Table [3J we notice that the 
and i?oo ratios are systematically below the data, while the absolute branching 
fractions for the colour- allowed modes with only charged particles in the final state are 
consistently above. A plausible interpretation of this trend is that A# is in fact smaller 
than the default value Xb = 350 MeV (a smaller value enhancing the colour-suppressed 
amplitude), and that the heavy-to-light form factors are about 10% smaller than the 
QCD sum rule central values that we assumed in Table [Hand used for "Theory I"@ 
Since we have already seen in Section [6.3.21 that a large value of As is in conflict with 
data, we adopt the "small Xb and form-factor hypothesis" and recalculate the theoretical 
prediction with the modified parameter values and ranges f+*(0) = 0.23±0.03, Aq P (0) = 
0.28 ± 0.03, A B (1 GeV) = (0.20±g;^) GeV. The result is column "Theory II" in Table EJ 
We consider these as our reference predictions to be tested against future more accurate 
experimental results. 

Comparing now "Theory II" with measurements, we find good agreement for the 
branching fractions, except for 7r°7r° as already discussed. However, the theoretical un- 
certainty is also large for this mode. The same feature is reflected in the ratios, where 
Rqq stands out as too small. The ratios Rq, R p c of colour-suppressed modes do not fit 
very well either, reflecting the fact that factorization likes to have [p°p°}l to be larger and 
7r°7T to be smaller than data. Nevertheless, uncertainties taken face value, there is no 
disagreement. The up ratios Ris are in good agreement with data, though the central 
value of i?5 shifts away from the data in scenario "Theory II" . Since R^ and R$ provide 

7 BrAv denotes the CP-averaged branching fraction. 

8 Alternatively, a smaller value of \V u b\ might be considered. Indeed, our central values imply 
\V u b\f+ "'(O) = 9.34 • 1CT 4 , which is slightly larger than J60]). 
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access to the real part of a<i (ftp) and ^(pn), respectively [5], precise measurements of 
the np final states should provide further insight into the mechanism that generates the 
colour-suppressed tree amplitude. 

7 Conclusion 

We computed the two-loop vertex corrections to the colour- allowed and colour-suppressed 
tree amplitude in QCD factorization, completing the calculation of these amplitudes at 
NNLO. Technically, the calculation amounts to a matching calculation from QCD to 
SCET, involving two-loop renormalization and infrared subtractions of evanescent oper- 
ators, and massive two-loop vertex integrals that depend on one dimensionless parame- 
ters, and the charm quark mass in some cases. We obtain fully analytic expressions for 
both amplitudes after integration over the Gegenbauer expansion, including the exact 
dependence on the charm quark mass. The massless result is in complete agreement 
with a recent independent calculation [13], as well as the charm mass dependence of the 
imaginary part [12], while the charm mass dependence of the real part agrees numerically. 

The NNLO vertex correction to the colour-suppressed tree amplitude is sizable, rang- 
ing from 10% to 25% for the real and imaginary part, respectively, and is a few percent 
for the colour-allowed amplitude. When combined with the 0(a 2 s ) correction to spectator 
scattering, already known from [6-8], we find a large cancellation, both in the real and 
imaginary parts, so that the overall NNLO correction to the topological tree amplitudes 
is small. This is somewhat unfortunate, since an enhancement rather than a cancellation 
in the colour-suppressed tree amplitude might have helped to cure the large discrepancy 
with experimental data in the channel. 

A dedicated phenomenological analysis of tree-dominated B decays to (quasi) two- 
body final states with pions and rho mesons shows that the QCD factorization approach 
obtains strong support from the factorization test performed in Section 16.3.11 Overall, 
the data is described very well within theoretical and experimental uncertainties, espe- 
cially for low values of A# ~ 200 MeV and smaller form factors. The most problematic 
observables remain those related to the 7r°7r° branching fraction, which is predicted too 
low. Since the combined NNLO correction is rather small, the generic features of factor- 
ization are unchanged compared to the NLO analysis of [5]. Numerical differences arise 
primarily from modified parameter choices. Our final results are the columns labelled 
"Theory II" contained in Tables [2] and [31 

Acknowledgement 

We would like to thank Bernd Jantzen for useful discussions, and Oleg Tarasov for 
collaboration at an early stage of this work. This work is supported in part by the 
DFG Sonderforschungsbereich/Transregio 9 "Computergestiitzte Theoretische Teilchen- 
physik". X.-Q. Li acknowledges support from the Alexander- von-Humboldt Stiftung. 
T.H. acknowledges support from the German Federal Ministry of Education and Re- 



31 




Figure 7: Selected master integrals. All momenta indicated are assumed to be incoming. 
Bold lines denote lines of mass m&, and thin lines are massless unless otherwise indicated. 
Zig-zag lines have virtuality um^. Dots on lines denote squared propagators, and the 
diagrams in question stand for both the undotted and the dotted integral, n represents 
a generalized propagator power. 
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A Master integrals 



In this Appendix we list a few master integrals on which our results add to the ones given 
in [11], for instance by giving closed forms valid to all orders in e = (4 — D)/2. They 
are depicted in Figure [7J The expansions in e are conveniently done with the package 
HypExp [36,37] and agree with the results given in [11]. Our notation for the integration 
measure is 

' d D k 



[''*]= l^-n- (68) 
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and we define the prefactor 



(69) 



(47r) D/2 r(l-e) ' 

All momenta are incoming, and the kinematics is such that q+qz+q^ = with q 2 = q 2 = 
and q\ = m\. We tacitly assume that all propagators in the integrals below contain an 
infinitesimal +ir). We then have 
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(71) 



For the next integral it is convenient to split up the integration over y into the three sub- 
intervals . . . u/(u + 1) . . . u . . . 1. Alternatively, one can compute this master integral 
by using the method of differential equations. 

1 



'5,3 



J [dh] J [dk 2 ] 



[{h - <? 4 ) 2 - m 2 ] (fci + k 2 + uq) 2 {h - g 3 ) 2 k\ k\ 



S 2 (m 2 b ) 



-i-2e r 2 (i-6)r(i + 26)r(i + 6)r(- e ) 

2e 2 



dy 



y 



-Ac 







u-y 



x 



2Fl(£ , 2£ ; ! - e ; V(u + ir,-v) ) _ 2e ; 1 — e ; ^t±M ) 



2 /_2\-l-2e 



-s£ K) 



-i ln 3 ( M ) + 1 ln(l - u) \n 2 {u) - \-n 2 ln(u) + Li 3 (l - u) 
o 2 o 
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fl n 2 

+Li 3 (w) + in ( - ln 2 (w) + Li 2 (l - u) + — 

+ - ^ ln(l - u) \n\u) - 4 ln 2 (l - u) ln 2 («) + ^Li 2 («) \n 2 {u) - -7t 2 ln 2 {u) 

8 6 2 3 

19 5 

+— vr 2 ln(l - u) ln(u) - 5 ln(l - u)Li 2 (u) ln(w) - HLi 3 (l - u) \a(u) - -Li 2 (w) 

O Zi 



-llLi 3 (u) ln(u) + Li 4 (l - u) + 13Li 4 (u) - 6S 2 , 2 (u) + 6 ln(l - u)C(3) - 

13 /5 
H 7r 2 Li 2 (w) - 6 ln(l - w)Li 3 (w) + ivr I - ln(l - u) ln 2 (w) - Li 2 (w) ln(tt) 



209tt 4 
360 



(72) 



h,i = J [dki] J [dk 2 ] J^TZ 
= S 2 (ml 



uq) 2 - m 2 ] ( kl + q 3 + uq) 2 (k 1 + k 2 ) 2 k\ 



2,2,-2e T 5 (l- e ) 



r(2 - 2e) 



x < — u 



r(3e - i)r(i - 2e) 

T(2 - 2e) 
r(l-3e)r(3e 



e 3?rie 2 Fi(l - e ,l-2e; 2 - 2e ; u) 



r 2 (l -e)r(2-3e) 



(73) 



{0,l-e,l-2e} , {} 
{0,l-2e} , {l-3e} 

The Meijer-G functions can be written as a linear combination of hypergeometric func- 
tions [97]. The latter are then expanded in e with HypExp. The next integral is a master 
only for n = 1 but for convenience we generalize one of the propagator powers. 

h,2 = J [dh] J [dk 2 ] 1 



[(h + g 4 ) 2 - m 2 ] n [(hi + k 2 + q A ) 2 - m 2 ] (h - q 3 ) 2 k 2 

S$ (-1)-" (m 2 ) J — — -OVA I 



^2 / i\-n I 2\l-n-2e F 2 (l / f > 

TV \ 44 

T(n) 



{1 — e — n, 1 — n, e, 2 — n — 2e} , {} 
{0, 1 — e — n} , {— n, e — 1} 



(74) 



/4 )3 = y [dfci] j [dk 2 \ 



i 



'4,4 



j\dHj [dk2] 



{(k 2 + q 4 ) 2 - m 2 } k\ (h + k 2 - qf (k 2 - uq) 2 

2 Fi(l,2e; 2-e; u) 
1 



_<* (rn 2,-2e r 4 (i-6)r( e )r(26)r(i-26) 

^ rlmJ r(2-2e)r(2-e) 



(75) 



[(ki + q 4 ) 2 - mf] (ki - uq) 2 \k\ - m 2 ] [(k x + k 2 ) - m 2 ] 
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1 1/5 uln(u)\ 2«Li 2 (l-«) 2uln 2 (u) 



1 " ^ 2e 2 e\2 u-1 J u-1 u-1 

5uln(u) it 2 19 „ . M/ . 1 . . 

-- ^rr + Y + Y + /44(^i-«) + o(6)j , (76) 

with 



/44(2, z) = - 



dwi f dw 2 



2m J 2m 

c-^—ioc C2 - 



— x)™ 2 ^ 1 csc(7ru;i) csc(7rw 2 )r(— 2^! — w 2 + l)T(wi + w 2 ) 
(w x - l)iuir (I - Wx) 

(77) 

and z = ml /ml. The two-fold Mellin-Barnes integration is along straight lines parallel 
to the imaginary axis, hence the real parts along the curves are constant. They read, 
respectively, c\ = 2/3 and c 2 = —1/2. The form (ITS"]) of enables us to obtain fully 
analytic results for the amplitudes a,\ and a 2 , since we can interchange the integration 
over u with the Mellin-Barnes integrations. 



h,5 = J [dh] J [dk 2 ] 



1 



[(h + q A ) 2 - m 2 ] (kx - uq) 2 \k\ - m 2 ] [(hi + k 2 ) - m 2 ] 



e (tt — 1) it — 1 u — 1 



(78) 



with 



c^+ioo C2 +? 00 



Ci —7 OO 



X 



27TZ J 27T2 

2 2mi - 1 7r 5 / 2 (l - x) w *z Wl csc(vru;i) csc(7rw 2 )r(-2wi - w 2 )r( Wl + w 2 + 1) 

io x r (| - u>i) 

(79) 

and Ci = —1/24 and c 2 = —5/6. 

J [dh] j [dk2] k\ (h + k 2 + q) 2 (k 2 - q 3 ) 2 (k 2 + uq) 2 

9 2 (m 2 ^ r 5 (l~e)r(e)r(2e)r(l-2e) 2 _ . . 

~^ K) r(2-2 e )r(2-3e)r(l + e) 6 2^1(1,26,1 + 6, «). (80) 
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i.e 



In the next integral we again generalize one of the propagator powers. 

1 



'4,7 



s 2 



[dki] J [dk 2 
i .-1) 



m 



[(k. + q^f ih-uqfikl-m^ [(h 

T2 l-2n-2e 



mi 



2\l-n-2e 



r 2 (i - e) 



r (n) 



u 



{1 - e - n, 1 - n, 0, 2 - n - 2e} , {} 



{0,1 



n 



Ml 



71 



e, e 



(81) 



The remaining integrals from Figure [7] can be found in the literature and shall not be 
given explicitly here: The integral l6,i through order O(e ) can be found in [98], whereas 
an all-order result for 1^2 is contained in [99]. The three-line master integral can be 
found in various papers [100-103]. 



B Auxiliary functions 
B.l Hard-scattering kernels 

In this Appendix we list the explicit expressions for the right-insertion hard-scattering 
kernels. The wrong-insertion kernels can be obtained from linear combinations of the 
right-insertion ones, see (1361) . The one-loop hard-scattering kernels in the CMM operator 
basis read 



n 



(i) 



Cf_ 

2N r 



-6L + ln 2 (l -u) + 2 ln(it) ln(l - it) - 
+ 31n(l 



21n(l - it) 



u 



u 



, 9/ , ln(it) , . , . . 7r 

ln 2 ( M ) + + 3 l n ( M ) + AU 2 {u) - — 

1 o 



22 



it 



+ in 



21n(l-u) -21n(it) -3 



T, 



(i) 



with 



L = ln 



mi 



(82) 



(83) 



At two loops, the kernels are conveniently subdivided into the following building 
blocks 



(2) 



1 (2),im 



1 (2),im 



(2),re , ■ rp(2),im\ 



,(2) 



T-i(2),re . ■ r7-i(2),«m 



(84) 
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Figure 8: Diagrammatic representation of the charm-mass dependent term. Three fur- 
ther diagram differences are not shown. 

The charm-quark contribution Tf (T^' re + iirT^' m ) to arises from charm-loop 
insertions into the gluon propagator of the non-factorizable one-loop diagrams. It is 
designed such as to give the difference between the contribution of a quark of mass 
m c and a massless quark, see Figure [HJ This construction ensures that we set ni = 4 
irrespective of whether we treat the charm quark as massive or massless. The explicit 
expressions for the two-loop hard-scattering kernels are 

T (2),re _ (47m 5 - 278u 4 + 1223m 3 - 2316m 2 + 2036u - 652) ln 4 (l - u) 
1 ~ 162(m- 1) 2 u 3 

(2m 3 + 4m 2 + 173m + 16) ln 3 (l - u) (4u 3 - 61u 2 - 436m + 16) ln 2 (l - u) 
81m 54m^ 

2 (73m 5 + 38m 4 - 1103m 3 + 2316u 2 - 2036u + 652) ln(u) ln 3 (l - u) 
+ 81(m- 1) 2 m 3 

(17m 3 + 300m 2 - 1098m + 978) \n 2 (u) ln 2 (l - u) 



7T 2 (9m 5 + 166m 4 - 1167m 3 + 2316m 2 - 2036m + 652) ln 2 (l - 

C1 In, VY273, 



27m 3 

~ ' 2036m + 652) In 2 (1 - u 

31(m - 1) 2 m 3 

(2M 5 -20M 3 + 125u 2 -76M-52)ln(M)ln 2 (l-M) 2, , 
+ ~ 27(,-l)% L + 9 ln <"> " "> 

7(u-2) 2 ln(2-ti)ln 2 (l -u) 16 T .,.,,„ 

+ - — 9(«-d» — 1 + y L,2(u) ln (1 - u) 

(2m 6 + 4m 5 - 191m 4 - 167m 3 + 1022m 2 - 646m - 6) \n 2 (u) ln(l - u) 
+ 27(m - 1)m 3 

tt 2 (2u 5 + 355m 3 - 623u 2 + 385m - 140) ln(l - u) 
81(m - 1) 2 m 

(4m 4 - 638m 3 + 1487m 2 - 1597m + 664) ln(u) ln(l - u) 
27(m - 1)m 2 
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U(u - 2) 2 Li 2 (w - 1) ln(l - u) 16 (6m 2 - 16m - 5) Li 3 (u) ln(l - u) 
+ 9(m- l) 2 + 27(m- l) 2 

2 (94m 3 - 271m 2 + 166m + 32) Li 2 («) ln(l - m) (1601m - 1172) ln(l - m) 



+ 



27(m - 1) 2 m 54m 
4 (4m 3 - 50m 2 + 183m - 163) ln(u)Li 2 (u) ln(l - u) 

(2m 3 - 436m 2 + 657m - 332) 1ii 2 (m) 8 (3m 2 - 14m - 19) C(3) ln(l - u) 



27(m-1)m 27(m-1) 2 

2 (20m 5 - 94m 4 + 292m 3 - 579m 2 + 509m - 163) L1 2 (m) 2 
27(m- 1) 2 m 3 

7T 2 (4m 4 - 435m 3 + 3174m 2 - 5346m + 2688) 64 T . 

+ — Li 3 (1-m) ln(l-M) 



162(m-1)m 2 9 

7T 4 (225m 5 + 378m 4 - 5914m 3 + 11598m 2 - 10198m + 3266) 
2430(m - 1) 2 m 3 

16 (m 3 + 2m - 2) HPL({-2, 2}, 1 - u) (2m 3 + 4m 2 - 313m + 273) 1ii 3 (m) 



27m 3 81(m-1) 

_ 8(m - 2) (5m 2 - 6m + 2) HPL({-1, 2}, 1 - u) 47 1ii 4 (m) 
27(m - 1)m 2 162 

tt 2 (m-2)(53m 3 -104m 2 -16m + 4)1ii(2-m) 65 2l 2/ , 
+ 81(m-1) 2 m 2 81" ln {U) 



7T 



2 (2m 5 + 4m 4 - 427m 3 + 1209m 2 - 2140m + 1320) 1ii(m) 



81(m- 1)m 2 

(1601m- 942) 1ii(m) 4 (m 2 - 3) ln(u)Li 2 (-u) 4vr 2 (m 3 + 2m - 2) Li 2 (M - 1) 



54(m - 1) 27m 3 81m 3 

2tt 2 (33m 3 + 400m 2 - 1472m + 1312) 1ii(1-m)1ii(m) 16 , f 
+ ^73 -HPL({-2,2},m) 

2tt 2 (2m 4 + m 3 + 5m 2 + 9m - 9) 1ii(m + 1) 16(m + 1) 2 HPL({-1, 2}, u) 
+ 81(m- 1)m 3 + 27(m- 1)m 

(8m 4 - 937m 3 + 1523m 2 - 1509m + 664) L1 2 (m) 2 (m 2 - 3) 1ii 2 (m) 1ii(m + 1) 



27(m - 1)m 2 27m 3 

2tt 2 (10m 5 - 208m 4 + 1183m 3 - 2340m 2 + 2060m - 660) U 2 (u) 

81(m- 1) 2 m 3 

2 (34m 4 + 530m 3 - 1551m 2 + 972m + 6) ln(M)Li 2 (M) 4 2 



27(m-1)m 3 81 
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+ — 7T Li 2 (-M) 



2 (4m 7 + 2m 6 - 39m 5 - 63m 4 + 678m 3 - 959m 2 + 362m - 6) Li 3 (l - u) 
+ 27(m- 1) 2 m 3 

16(6M 3 + 50M 2 -183M + 163)ln(M)Li 3 (l-M) 68 2 



+ — ln 2 (M)Li 2 (M) 



27m 3 27 

14(m - 2) 2 L1 3 (m - 1) 4(M 2 -3)Li 3 (-M) 104 

+ ^ — ln(u)Li 3 («) 



9(m - l) 2 27m 3 9 

2 (4m 7 + 2m 6 - 39m 5 + 743m 4 - 2747m 3 + 3350m 2 - 1292m - 6) L1 3 (m) 

27(m - 1) 2 m 3 

4 (86m 5 - 406m 4 + 1553m 3 - 2919m 2 + 2569m - 823) Li 4 (l - u) 

27(m - 1) 2 m 3 

4 (96m 5 - 372m 4 + 1283m 3 - 2316m 2 + 2036m - 652) Li 4 («) 

27(m - 1) 2 m 3 

16 (5m 5 + 36m 4 - 289m 3 + 579m 2 - 509m + 163) Li 4 (^) 

27(m - 1) 2 m 3 

2 (62m 5 - 397m 4 + 1012m 3 - 1112m 2 + 428m + 2) C(3) 
9(m- 1) 2 m 3 

8(3m 3 + 100m 2 - 366m + 326) 1ii(m)C(3) 16tt 2 (2m - 1) ln(2) 3793 



27m 3 27(m - l)u 36 

T f 32 n 3/, \ 2(27m + 4)1ii 2 (1-m) 32 , . . , 2 . 
+L (-- ln 3 (l - m) + J ±-± 1 - - ln(M) ln 2 (l - u) 

32, 2 , Nl „ , 2(13u - 10) ln(l - u) 4(51u - 4) \n(u) ln(l - u) 
1ii 2 (m) ln(l -u) + — — + — — — — 

27 V ; V ; 3m 27m 

64 T . . Nl M . 32 2l M , 321ii 3 (m) 2(51m - 43) 1ii 2 (m) 

Li 2 (m) ln(l - u) H 7T 2 ln(l - u) H — ; '—— — 

27 K J K J 81 V ; 81 27(m-1) 

2tt 2 (27m - 23) 26(3m - 2) 1ii(m) 8 (39m 2 - 39m + 2) Li 2 (u) 



81(m-1) 9(m-1) 27(m-1)m 

64 . \ t • / \ 64 t- /, \ 64Li 3 (M) \ 26L 2 
-- In(«)Li 2 («) - -Li 3 (l - u) + 2 y > - 55 J - — , (85) 

where the functions HPL denote the harmonic polylogarithms [43]. We proceed with 

T (2),im 8 3 2 (M 3 + 2M 2 -53M + 10)ln 2 (l-M) 8 2 

1 = 81 ^ ^ 27m" + -ln(M)ln (1 -m) 

(4m 2 + 455m - 40) ln(l - u) 8 . . . . 64 2 
+^ " + 27 W«) ln(l - u) - ln(l - m) 
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2 (2m 4 + 4m 3 - 40m 2 + 31m - 12) ln(u) ln(l - u) 56, , H 
~- 27(m-1)m " " 27 ln {U) ln(1 " M) 

2tt 2 (3m 3 + 5m 2 - 8m - 5) (2m 3 + 4m 2 + 38m - 57) In 2 (m) 4 In 3 (m) 
+ 81(m- 1) + 27(m - 1) + 27 

(4m 2 + 221m-210)1ii(m) 40 2l . . 4 T . . . 8 , . . T . . . 

+ ^tTT M«) - — mLi 2 (m) - - 1ii(m)Li 2 (m) 



27(m- 1) 81 v ' 27 v ' 3 

20Li 2 (M) 8Li 2 (M) 4L1 2 (m) 80 t . , 40Li 3 (M) 1609 
+ 277m^-^7m- + ^^-27 L13(1 - m) + ^ 5T 

/_ 16 2 4(39m-35)1ii(m) _ 64Li 2 (M) 16tt 2 

+ V 27 11 _ 27(m - 1) 27 +_ 8T 

4(39m - 4) ln(l - u) 32 , . . , .„ . 161ii 2 (m) 26 \ 

+ J ^ i - 27 K«) Ml - «) + - y • (86) 



ln = (u) ln(1 _ u) _ - 32) tad - «) 



27 v ' y 81m 

8(19m - 16) 1ii(m) ln(l - u) + 16^ 2 ^ + 8(14m - 11) 1ii 2 (m) 



81(m-1) 27 v ' 81(m-1) 

4tt 2 (19m-16) 4(51m - 32) 1ii(m) 16 81n3 ( M ) 



243(m- 1) 81(m - 1) 27 v ' 27 

8Li 2 (M) 16Li 2 (M) 232Li 2 (M) 16 . 16Li 3 (w) 250 
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27(m - 1) 27m 81 27 ° v ' 27 27 

r ( 8 , 2m n 8(3m - 2) ln(l - m) 16, . Nl „ , 
+L ("27 ln (1 " U) ~ 27* " " 27 HU) ln(1 " U) 

81ii 2 (m) 8(3m-2)1ii(m) 32L1 2 (m) 8tt 2 136 \ 8L 2 
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^ = l Ma(1 _ tt) _lM Ml _ u) _^ + 801^ + 68 
i,n, 27 V ; 81 V ; 27 81 27 

-^L ^21n(l -u) -21n(u) -3^ . (88) 
The terms arising from the bottom quark-loop insertion into the gluon line read: 
T?r = L (-1 ln 2 (l - u) - 8 ^-2) ln(l-M) _ 16 _ 

1,7 > V 27 27 « 27 
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with the abbreviations 



u 1—u 

"1 = irtr , 0i = t^TT ' (91) 



and the auxiliary functions 

Pl (x) = 12 In 2 ( y^±L ) - 3 ln 2 (l - x) + 12 In ( H^±}1) ln (i 
V2vl — V v 1 — £ / 

-24 In (^=) In(^) - 24Li 2 - 24Li 2 (-^) 
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* w = 21n tfn^J ni7r=^J - 3L HV5TI 1+21n (2) ' 



p 3 (x) = In 



The charm-mass dependent parts mentioned at the beginning of the section read 
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where now 



M 1 — M 

4z j. — 4z 



= ^tt ' *' = nfrr ' (95) 

4z "r 1 4z ^ 1 

as well as 2 = m 2 /m 2 . The other functions, notably / 44 , / 45 , g 7 , and g 8 , can be found in 

_____ lcy\ 

(JTTj) . (1791 . and (199]) . respectively. The two-loop expressions for the kernel T 2 related to 
the insertion of the operator Q 2 are simpler. They are 
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2 (m 3 + 2m 2 + 7u - 4) ln 3 (l - m) (2m 3 - 27u 2 - I2u + 8) ln 2 (l - u) 



+ 



27m 9m 2 
16 (2m 5 - 8m 4 + 16m 3 - 36m 2 + 35m - 12) ]n(u) ln 3 (l - u) 41n 4 (M 



27(u - l) 2 u 3 27 
4 (u 3 - 6m 2 + 33m - 36) ln 2 ( M ) ln 2 (l - u) 2{u - 2) 2 ln(2 - u) ln 2 (l - u) 



9m 3 3(m - If 



in 



2 f u 5 + 2u 4 - 13m 3 + 36u 2 - 35m + 12) ln 2 (l - u) 



27(m - 1) 2 m 3 9 



- -U 2 (u) ln 2 (l - u 
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2(w 5 -23w 3 + 87« 2 -92w + 30)ln(w)ln 2 (l-w) 8 , „, „ 

; tt; In M m l — u) 

9(m-1) 2 m 27 V ; V ; 

2(m 6 + 2m 5 -5m 4 + 12m 3 -63m 2 + 51m-3)1ii 2 (m)1ii(1-m) 4 2 2 



9(m - 1)m 3 9 
(151m - 86) ln(l - u) 2tt 2 (m 5 - 8m 3 + 75m 2 - 95m + 30) ln(l - u) 



-7r 2 ln 2 (M) 



9m 27(m - 1) 2 m 

32tt 2 (m 3 + 2m 2 - 11m + 12) 1ii(m) ln(l - u) 4(m - 2) 2 Li 2 (M - 1) ln(l - u) 



+ 



27m 3 3(m - l) 2 

2 (2m 4 - 50m 3 + 14m 2 + 81m - 42) 1ii(m) ln(l - m) 16 T . „ , , „ 

Li 3 (l — u) m(l — u) 



9(m-1)m 2 3 
4 (12m 3 - 53m 2 + 58m - 20) Li 2 (u) ln(l - u) 32(2m + 1)C(3) ln(l - u) 



9(m-1) 2 m 9(m-1) 2 
8 (4m 3 - 2m 2 + 11m - 12) ln(M)Li 2 (M) ln(l - u) 2 (m 2 - 3) 1ii 2 (m) 1ii(m + 1) 



9m 3 9m 3 
4 (8m 5 - 16m 4 + 24m 3 - 36m 2 + 35m - 12) L1 2 (m) 2 2tt 2 (m 2 - 3) 1ii(m + 1) 



9(m - 1) 2 m 3 9m 3 
tt 2 (2m 4 + 3m 3 - 100m 2 + 298m - 168) 2 (m 3 - 34m 2 + 7m + 21) 1ii 2 (m) 



27(m-1)m 2 9(m-1)m 
tt 4 (62m 5 + 2m 4 - 195m 3 + 720m 2 - 700m + 240) 4 (m 2 - 3) ln(u)Li 2 (-«) 



405(m - 1) 2 m 3 9m 3 

27t 2 (m-2) 2 1ii(2-m) _ (151m- 102) 1ii(m) 2(m - 2) (m 2 + 4m - 3) 1ii 3 (m) 
+ 3(m- l) 2 9(m- 1) + 27(m - 1) 

2vr 2 (m 5 + 2m 4 - 23m 3 - 15m 2 + 136m - 96) \n(u) 4(m - 2) 2 Li 3 (M - 1) 



27(m-1)m 2 3(m-1) 2 

4tt 2 (4m 5 - 22m 4 + 61m 3 - 144m 2 + 140m - 48) Li 2 (u) 
27(m- 1) 2 m 3 

4(6M 4 -26M 3 + 97M 2 -75M + 3)ln(M)Li 2 (M) 8 2 

g(u _ i)u3 + - In (M)Ll 2 (M) 

4 (2m 7 + m 6 - 18m 5 + 79m 4 - 111m 3 + 63m 2 - 10m - 3) Li 3 (l - u) 

9(m- 1) 2 m 3 

32(m - 4)(2m - 3) ln(M)Li 3 (l - u) 2 (4m 4 - 97m 3 + 51m 2 + 83m - 42) Li 2 (u) 



9m 3 9(m - 1)m 2 
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_4 (2m 7 + m 6 - 12m 5 - 12m 4 + 139m 3 - 220m 2 + 102m - 3) Lis(u) 

9(m - 1) 2 M 3 

16, , NT , N 8(2M 5 + 14M 4 -69M 3 + 180M 2 -175M + 60)Li 4 (l-M) 

-T^'W") + " ^T)^ — 1 

8 (2m 5 - 22m 4 + 57m 3 - 144m 2 + 140m - 48) U A (u) 1507 
+ 9(m - 1) 2 m 3 + ^8~ 

32 (4m 4 - 14m 3 + 36m 2 - 35m + 12) Li 4 (^) 32(m - 4) (2m - 3) 1ii(m)C(3) 



9(m - 1) 2 m 3 9m 3 

4 (36m 5 - 86m 4 - 67m 3 + 227m 2 - 110m + 3) C(3) 4 (m 2 - 3) Li 3 (-«) 
9(m - 1) 2 m 3 <hfi 

r ( 4 t 2/, x 4(3m - 2) ln(l - m) 8. , N . M 
+L [-- ln 2 (l - u) - -± '- - - \n{u) ln(l - m) 

+ 41^(m) _ 4(3m - 2) ln(M) _ Ig^tt) + 1^ + 34) + 4L 2 , (96) 

if- = — ln 3 (l - m) - 2(M 3 + 2M 2 -5M + 6)ln 2 (l-M) 16 _ 
2 27 v ; 9m 9 y j y j 

4 (m 2 + 29m - 5) ln(l - m) 4 (m 4 + 2m 3 - 5m 2 + 6m - 2) 1ii(m) ln(l - u) 



9m 9(m - l)u 

+ ^fLi 2 (M) ln(l - u) - V ln(l - u) - ^ 1ii 2 (m) ln(l - u) 
9 9 9 

+ 81ii 3 (m) 2(m 3 + 2m 2 + m- 1)1ii 2 (m) 2tt 2 (3m 3 + 5m 2 - 11m + 4) 



27 9(m-1) 27(m-1) 

4, s n , s 8 9l . , 4 . . 32 / \ T / \ 4Li 2 (M) 
+ -(m + 21) 1ii(m) - — vr 2 1ii(m) + -mLi 2 (m) - — 1ii(m)Li 2 (m) + 



9 V ' v y 27 v y 9 zv y 9 v y zv y 9(u-l) 

8Li 2 (M) 4Li 2 (M) 16 T . , 16Lia(u) 155 
- ; +— Li 3 (1-m) + 7r~^ + 



9m 3 9 JV ' 9 9 

+L (-| M l-„) + iM^ + 4). (97) 

B.2 Charming functions 

In this Appendix we list the expressions for the functions that depend on m c . We define 
z = ml/ml and 

Az 

V= (v/i^TT+lf 

45 



K = 1 + \/z . 



(98) 



In terms of these variables our building blocks read 



qi (z) = ln 3 ^) + 6 111(77) Liafa) - 6 14(77) + 6((3) , 
q 2 (z) = In^V^z + 1 , 

q 3 (z) = 9 ln 4 (r7) - 12tt 2 In 2 (77) + 144 Li 2 (7?) - 48tt 2 L1 2 (t7) + 4vr 4 , 

q 4 (z) = In 3 (77) - 27T 2 ln(77) + 12Li 3 (?7) - 12C(3) , 

q 5 (z) = - ln 3 (77) - 4 ln 2 ^) 14(77) + 2tt 2 111(77) " 12 U 3 ( V ) + 12((3) , 

q 6 (z) = [-3 ln 2 (77) - 12 Li 2 (r7) + 2tt 2 ] V^+T , 

q 7 (z) =ln 2 (2)+2Li 2 (l-^) + 7r 2 , 

q 8 (z) = 2 Iu(k) ln(z) - ln(l - z) ln(z) - 4 Li 2 (V*) + U 2 (z) + tt 2 , 

q 9 (z) = Li 3 (l - 2) - 214(1 - v^) - 2Li 3 ( v / i) + Li 3 (z) +2^3^^^-^) - 2((3) , 

q w (z) = - ln 3 (z) + 3 ln(l - z) \n 2 (z) + 6 ln 2 (/t) ln(z) - 4tt 2 \n(z) - 4 ln 3 (/t) 
+87r 2 ln(/t) + 12g 9 (z) , 

q n (z) = ln\z) +241n 2 (/t)ln 2 (z) -481n(l - 2) ln(/t) ln 2 (^) +87r 2 ln 2 (^) 

+ 12 Li 2 (z) ln 2 (z) - 8tt 2 ln(l - z) ln(z) + 16tt 2 1ii(k) \n(z) + 16 ln 3 (/t) ln(z) 
-96 ln(l - y/z) \n{z) Li 2 (-y/z) - 192 Li 2 (-^) Li 2 (V*) + 40tt 2 Li 2 (z) 
-96 ln(re) ln(z) Li 2 (^i) - 96tt 2 Li 2 (^i) + 8tt 4 - 48 ln(z) g 9 (z) , 

q 12 ( z ) = -4 ln 3 (l - z) + 6 ln(z) ln 2 (l - 2) + 12 ln(2/t) ln 2 (l - z) + 3 ln 2 (^) ln(l - z) 
-12 ln 2 (2/t) ln(l - z) - 24 ln(z) ln(/c) ln(l - 2) - 4tt 2 ln(l - z) + 4 ln 3 (/t) 
+ 18 ln 2 (/t) ln(z) - 12 ln 2 (2) ln(z) + 8tt 2 1h(8*/k) - 6 ln(/t) ln 2 (^) 

+24 ln(2) \n(2z) ln(re) - 12 ln(z) Li 2 (z) - 24 ln(z) Li 2 {^j - 24 Li 3 ( - 




+24 Li 3 (1 - yfz) - 72 Li 3 (y/zj + 48 ln(z) Li 2 (y/z) + 12 Li 3 (z) - 24((3) . (99) 




+24Li 3 
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